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We investigate the vacuum expectation value of the surface energy-momentum tensor for a massive 
scalar field with general curvature coupling parameter obeying the Robin boundary conditions on 
two codimension one parallel branes in a (D + l)-dimensional background spacetime AdSz>i+i X £ 
with a warped internal space E. These vacuum densities correspond to a gravitational source of the 
cosmological constant type for both subspaces of the branes. Using the generalized zeta function 
technique in combination with contour integral representations, the surface energies on the branes 
are presented in the form of the sum of single brane and second brane induced parts. For the 
geometry of a single brane both regions, on the left and on the right of the brane, are considered. 
At the physical point the corresponding zeta functions contain pole and finite contributions. For an 
infinitely thin brane taking these regions together, in odd spatial dimensions the pole parts cancel 
and the total zeta function is finite. The renormalization procedure for the surface energies and 
the structure of the corresponding counterterms are discussed. The parts in the surface densities 
generated by the presence of the second brane are finite for all nonzero values of the interbrane 
separation and are investigated in various asymptotic regions of the parameters. In particular, 
it is shown that for large distances between the branes the induced surface densities give rise to 
an exponentially suppressed cosmological constant on the brane. The total energy of the vacuum 
including the bulk and boundary contributions is evaluated by the zeta function technique and the 
energy balance between separate parts is discussed. 

PACS numbers: 04.62,+v, 03.70.+k, ll.10.Kk 



I. INTRODUCTION 

The braneworld scenario provides an interesting al- 
ternative to the standard Kaluza-Klein compactifica- 
tion of the extra dimensions. The simplest phenomeno- 
logical models describing such a scenario are the five- 
dimensional Randall-Sundrum type braneworld models 
(for a review see 0,0)- From the point of view of embed- 
ding these models into a more fundamental theory, such 
as string/M-theory, one may expect that a more complete 
version of the scenario must admit the presence of addi- 
tional extra dimensions compactified on an internal man- 
ifold. From a phenomenological point of view, the consid- 
eration of more general spacetimes offer a richer geomet- 
rical structure and may provide interesting extensions of 
the Randall-Sundrum mechanism for the geometric ori- 
gin of the hierarchy. More extra dimensions also relax 
the fine-tunings of the fundamental parameters. These 
models can provide a framework in the context of which 
the stabilization of the radion field naturally takes place. 
In addition, a richer topological structure of the field con- 
figuration in transverse space provides the possibility of 
more realistic spectrum of chiral fermions localized on 
the brane. Several variants of the Randall-Sundrum sce- 
nario involving cosmic strings and other global defects 
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of various codimensions have been investigated in higher 
dimensions (see, for instance, Q and references therein). 

Motivated by the problems of the radion stabilization 
and the generation of cosmological constant, the role of 
quantum effects in braneworlds has attracted great deal 
of attention A class of higher dimensional mod- 

els with the topology AdS^+i x S, where £ is a one- 
parameter compact manifold, and with two branes of 
codimension one located at the orbifold fixed points, is 
considered in Refs. H^,|23. In both cases of the warped 
and unwarped internal manifold, the quantum effective 
potential induced by bulk scalar fields is evaluated and it 
has been shown that this potential can stabilize the hier- 
archy between the Planck and electroweak scales without 
fine tuning. In addition to the effective potential, the in- 
vestigation of local physical characteristics in these mod- 
els is of considerable interest. Local quantities contain 
more information on the vacuum fluctuations than the 
global ones and play an important role in modelling a self- 
consistent dynamics involving the gravitational field. In 
the previous papers 0, [II we have studied the Wight- 
man function, the vacuum expectation values of the field 
square and the energy-momentum tensor for a scalar field 
with an arbitrary curvature coupling parameter obey- 
ing Robin boundary conditions on two codimension one 
parallel branes embedded in the background spacetime 
AdSd+i x £ with a warped internal space S. For an 
arbitrary internal space E, the application of the gener- 
alized Abel-Plana formula 1501 allowed us to extract form 
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the vacuum expectation values the part due to the bulk 
without branes and to present the brane induced parts in 
terms of exponentially convergent integrals for the points 
away from the branes. 

The braneworld corresponds to a manifold with bound- 
aries and the physical quantities, in general, receive both 
volume and surface contributions. In particular, the 
contributions located on the visible brane are of spe- 
cial interest as they are direct observables in the theory. 
In Ref. |35| the vacuum expectation value of the sur- 
face energy-momentum tensor is evaluated for a massive 
scalar field subject to Robin boundary conditions on two 
parallel branes in (D + l)-dimensional AdS bulk (for sur- 
face densities induced by quantum fluctuations of a bulk 
scalar field on a single brane in Rindler-likc spacctimes 
see ^^)- ^ has been shown that for large distances be- 
tween the branes the induced surface densities give rise 
to an exponentially suppressed cosmological constant on 
the brane. In the Randall-Sundrum braneworld model, 
for the interbrane distances solving the hierarchy problem 
between the gravitational and electroweak mass scales, 
the cosmological constant generated on the visible brane 
is of the right order of magnitude with the value sug- 
gested by the cosmological observations. In the present 
paper we investigate similar issues within the framework 
of higher dimensional braneworld models with warped 
internal spaces. 

The paper is organized as follows. In the next section 
we describe the corresponding bulk geometry and outline 
the structure of the eigenmodes and the Kaluza-Klein 
spectrum. An integral representation is constructed for 
the zeta function related to the vacuum expectation value 
of the surface energy-momentum tensor. This function 
is presented in the form of the sum of single brane and 
second brane induced parts. The analytic continuation 
of the zeta function for the model with a single brane 
is done in Section IlIII The surface energy- momentum is 
investigated on both sides of the branes and the renor- 
malization procedure is discussed. In Section llVl we study 
the surface energy density induced by the presence of the 
second brane. Various limiting cases are considered and 
the cosmological constant induced on the brane is esti- 
mated. The Section[V]is devoted to the evaluation of the 
total vacuum energy in the two-brane setup on the base 
of the zeta function technique. The relation between the 
bulk and surface energies is discussed. The last section 
contains a summary of the work. 



II. SURFACE ENERGY-MOMENTUM TENSOR 
AND THE RELATED ZETA FUNCTION 

Consider a massive scalar field (p(x) with curvature 
coupling parameter £. The corresponding field equation 
has the form (we adopt the conventions of Ref. [5l| for 
the metric signature and the curvature tensor) 

(V M V M + m 2 + (R) <p(x) = 0, (1) 



where Vm is the covariant derivative operator and R is 
the scalar curvature for a (D+l)-dimensional background 
spacetime. For the most important special cases of mini- 
mally and conformally coupled scalars one has £ = and 
C = Cd = (D — 1)/4D, respectively. We will assume that 
the background spacetime has the topology AdS£> 1+ i x £ 
and is described by the line element 

ds 2 = gMN<lx M dx N 

= e~ 2k °y (^dx^dx" - ltk dX l dX k ) - dy 2 , (2) 

where T}^ is the metric tensor for Z?i-dimensional 
Minkowski spacetime, kjj is the inverse AdS radius, 
and the coordinates X 1 , i — 1,...,Z?2, cover the in- 
ternal manifold S, D = D% + Z?2- The Ricci scalar 
corresponding to line element @ is given by formula 
R = -D(D + l)k 2 D - e 2kDV R h) , with R h) being the 
scalar curvature for the metric tensor 7^. In the discus- 
sion below, in addition to the radial coordinate y we will 
also use the coordinate z = e k ° y /kp, in terms of which 
the line element is written in the form conformally related 
to the metric in the direct product spacetime i?^ 1,1 ) x £ 
by the conformal factor (kjjz)^ 2 . 

We are interested in one-loop vacuum effects induced 
by quantum fluctuations of the bulk field f(x) on two 
parallel branes of codimension one, located at y = a and 
y = b, a < b. We assume that on the branes the field 
obeys Robin boundary conditions 

[Ay + Bydy^j lf(x) = 0, ?/ = <Z , 6 , (3) 

with constant coefficients A y , B y . This type of condi- 
tions is an extension of Dirichlet and Neumann bound- 
ary conditions and appears in a variety of situations, in- 
cluding the considerations of vacuum effects for a con- 
fined charged scalar field in external fields, spinor and 
gauge field theories, quantum gravity and supergravity. 
Robin boundary conditions naturally arise bulk fields in 
braneworld models. 

The branes divide the space into three regions corre- 
sponding to — oo < y < a, a < y < b, and b < y < oo. In 
general, the coefficients in the boundary conditions © 
can be different for separate regions. In a higher dimen- 
sional generalization of the Randall-Sundrum braneworld 
model the coordinate y is compactified on an orbifold 
S 1 /Z2, of ltngth I, —I ^ y ^ I, and the orbifold fixed 
points y = and y = I are the locations of two branes. 
The corresponding line-element has the form J2J) with the 
warp factor e _2fcD ' a '. In these models the region between 
the branes is employed only. For an untwisted bulk scalar 
with brane mass terms c a and c&, the corresponding ratio 
of the coefficients in the boundary condition ® is deter- 
mined by the expression (see, e.g., Refs. ll6L I38L for 
the case of the bulk AdS^+i and Refs. |26l l43l| for the 
geometry under consideration) 

B 3 2 w 
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In the supersymmetric version of the model |52( one has 
q, = — c a and the boundary conditions are the same for 
both branes. For a twisted scalar, Dirichlet boundary 
conditions are obtained on both branes. 

For the geometry of two parallel branes in AdSd+i x £ 
with boundary conditions @, the Wightman function 
and the vacuum expectation values (VEVs) of the field 
square and the bulk energy-momentum tensor are inves- 
tigated in Refs. O n manifolds with boundaries 
the energy-momentum tensor in addition to the bulk part 
contains a contribution located on the boundary. For an 
arbitrary smooth boundary dM with the inward-pointing 
unit normal vector n L , the surface part of the energy- 
momentum tensor for a scalar field is given by the for- 
mula p53| 

T$ N =6(x;dM) T MN, (5) 

where the "one-sided" delta- function 5(x;dM) locates 
this tensor on dM and 

tmn = C,ip 2 K M N - (2C - l/2)h M N<pn L V L ip. (6) 

In Eq. ©, Iimn = 9mn + ^m^at is the induced metric 
on the boundary and Kmn — h^jh^V ^np is the corre- 
sponding extrinsic curvature tensor. 

By expanding the field operator over a complete set 
of eigenfunctions {(p a {x), tp* a (x)} obeying the boundary 
conditions and using the standard commutation rules, for 
the VEV of the operator tmn one finds 

(0\t M n\0) = ^2T M N{<Pa(x),ip* a (x)} 7 (7) 

a 

where the bilinear form tmnI^P, "0} on the right is deter- 
mined by the classical energy-momentum tensor (0 and 
the collective index a can contain both discrete and con- 
tinuous parts. Below in this section, we will consider the 
region between the branes. The corresponding quantities 
for the regions y ^ a and y ^ b are obtained as limit- 
ing cases and are investigated in the next section. As 
we have mentioned before, in the orbifolded version of 
the model, the only bulk is the one between the branes. 
In this region the inward-pointing normal to the brane at 
y — j , j — a, 6, and the corresponding extrinsic curvature 
tensor are given by the relations 

n« M = n®5%, K% = -n^k D h MN , (8) 

where is defined in formula By using relations 
(JSJ) and the boundary conditions, the VEV of the sur- 
face energy-momentum tensor on the brane at y = j is 
presented in the form 

(0|rj$v|0) = -h MN n^k D C 3 (0\ V 2 \0)y= 3 . (9) 
with the notation 



From the point of view of physics on the brane, Eq. 
© corresponds to the gravitational source of the cos- 
mological constant type with the surface energy density 
= {0\t^°\0} (surface energy per unit physical vol- 
ume on the brane at y = j or brane tension), stress 
Pj — — (0|ti 1 |0), and the equation of state 



It is remarkable that this relation takes place for both 
subspaces on the brane. It can be seen that this re- 
sult is valid also for the general metric g^ a instead of 
77^0- in line element ijSjl. For an untwisted bulk scalar 
in the higher dimensional generalization of the Randall- 
Sundrum braneworld based on the bulk AdSot+i x 2 the 
coefficient in Eq. @ is given by the formula 

Oj = 4DC(C - Cd) + (12) 

In particular, the corresponding surface energy vanishes 
for minimally and conformally coupled scalar fields with 
zero brane mass terms. 

In order to evaluate the expectation value Q we need 
the corresponding eigenfunctions. These functions can 
be taken in the decomposed form 

= -% We (13) 
v /2a;^ > „(27r)- L ' 1 1 

where 

k^ = (wp.n.k), up,„ = ^k 2 +m 2 0n , k = |k|, (14) 

and rrifj^n are separation constants. The modes ipp(X) 
are eigenfunctions for the internal subspace: 

[A (7) + Ci? (7) ] MX) = -A^W, (15) 

with the eigenvalues Aj| and the orthonormalization con- 
dition 

j d D -X Vt^jPW ( x ) = 5 w ■ (is) 

In Eq. (|15fl . A/ 7 ) is the Laplace-Beltrami operator for 
the metric 7y. In the consideration below we will as- 
sume that A/3 ^ 0. The dependence on f3 and n in the 
separation constants is factorized: 

m l,n= m l + X % (17) 

where the mass spectrum m n is determined by the radial 
equation together with the boundary conditions. From 
the field equation (Q, one obtains the equation for the 
function f n (y) with the solution 



C i = C-(2C-l/2)I,/(fc I3J B J -)- 



(10) 



Mv) = C n e Dk °y^g^(m n z a ,m n z), 



(18) 
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where the normalization coefficient C n is determined be- 
low, and 

gi j \u,v) = J v (v)yW(u) - Y v (v)J<P(u), j = a,b. (19) 

Here J„ (x) , Y v (x) are the Bessel and Neumann functions 
of the order 



v = Jd 2 /4,- D{D + 1)C + m 2 /k%, (20) 

and the ^-coordinates of the branes are denoted by Zj — 
e DJ /hp, j = a,b. In formula l|19l) and in what follows 
for a given function F{x) we use the notation 

F^\x) =A j F(x) + B j xF'(x), j = a,b, (21) 

with the coefficients 

Aj = Aj + BjknD/2, Bj = Bjk D . (22) 

In the discussion below we will assume values of the 
curvature coupling parameter for which v is real. For 
imaginary v the ground state becomes unstable |H3.l55| . 
Note that for conformally and minimally coupled mass- 
less scalar fields one has v = 1/2 and v = D/2, respec- 
tively. 

The function l|18|) satisfies the boundary condition on 
the brane at y = a. Imposing the boundary condition 
on the brane y = b we find that the eigenvalues m n are 
solutions to the equation 

gi ab \m n z a ,m n z b ) = 0. (23) 

where the function on the left is defined by 

g ( « b \u, v) = Ji a \u)Y^{v) - Y^(u)Jl b \v). (24) 

Equation l|23[) determines the Kaluza-Klein (KK) spec- 
trum along the transverse dimension. We denote by 
u = 7^,1, n = 1,2,..., the positive zeros of the func- 
tion gi ab \u,uz b / z a ), arranged in the ascending order, 
7^,n < 7i/,7i+l- The eigenvalues for m„ are related to 
these zeros as 



m n = fcr>7^n e 



-kr>a 



^iv,ril Za • 



(25) 



From the orthonormality condition of the radial 
functions, 



dye^ D ^yf n (y)f nl (y) = <W, (26) 
for the coefficient C n in Eq. i|18|) on e finds 



7rrn n Y^ b \m n z b )/Y^ a '(m n z a ) 
k ° ^9 ( ^ b \uz a ,uz b )\ u=mn 



(a). 



el 



(27) 



Note that, as we consider the quantization in the region 
between the branes, a ^ y ^ b, the modes defined by i|18|) 
are normalizable for all real values of v from Eq. (|20|l . 



As it follows from formula J^J, the VEV of the surface 
energy-momentum tensor can be obtained from the VEV 
of the field square evaluated on the branes. The VEV 
of the field square on the bulk is derived from the cor- 
responding Wightman function in the coincidence limit. 
By using the generalized Abel-Plana formula for the 
summation over the zeros 7„. n , in Ref. ^3 this VEV is 
presented in the form 

(0|/|0) = <^><°> + (^)"' - 2kS- 1 z D l3o, 



J *0 



xJ2\MX)\ 2 duu(u* - \}) D ^-i 



xflj^(uz a , uz b )G^ 2 (uzj, uz) 



(28) 



where j — a and j — b provide two equivalent represen- 
tations and 



Pd 1 
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(4 7 r) D i/ 2 r(D 1 /2)" 
Here and in what follows we use the notations 



Q, av (u,v) 



Ki b \v) 



Ki a \u)Gr\u,v) 

li a \u) 
ll b) (v)Gl ab \u,vy 



(ab). 



(29) 



(30) 



(31) 



with the modified Bessel functions I u (u), K v (u), and 

G&(u,v) = I v {v)Ky\u)-K v {v)iy\u), (32) 
Gl ab \u,v) = Ki a \u)li b \v)-li a \u)Ki b Hv). (33) 

The first term on the right of Eq. J2HI, (^ 2 )^°\ corre- 
sponds to the VEV in AdSz^+i X X spacetime without 
branes, (<p 2 )^ is induced by a single brane located at 
y = j when the second brane is absent, and the last term 
is induced by the presence of the second brane. After 
the standard renormalization of the boundary free part, 
the VEV given by Eq. is finite in the bulk and di- 
verges on the branes. As a result, the VEVs of the field 
square on the branes cannot be obtained directly from 
this formula and an additional renormalization proce- 
dure is necessary. One possibility is to remove from l|28|) 
the terms which diverge for the points on the branes. 
An alternative way is to restrict from the begining to 
the fluctuations on the brane and to use the dimensional 
rcgularization or the generalized zeta function technique 
(for a discussion of the relation between different meth- 
ods see, for instance, Ref. ^3). Here we will follow the 
last method. 

Substituting the eigenfunctions (|13|l into the corre- 
sponding mode sum and integrating over the angular part 
of the vector k, for the expectation value of the energy 
density on the brane at y = j we obtain 



.00 



Jo 



dk k 



\ip0(X)\ 2 m n gu (m n zi,m„Zj 



yy 

n =l LOf3,n g^gl (UZ a ,UZ b ) 



■,(34) 
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where j, I = a, 6, and I ^ j, ujp^ n is denned by formulas 

(THjl . fTTjl . and we have used the relation <7^(u,w) — 
2Bj/n. To regularize the divergent expression on the 
right of this formula we define the function 



D Bj(3 Dl - 



3 



a JO 



m n gi l \m n zi,m n Zj) 

d (ab) i \i 

— g„ ■(uz ai uz b )\ 



(35) 



where an arbitrary mass scale p is introduced to keep the 
dimension of the expression. After the evaluation of the 
integral over fc, this expression can be presented in the 
form 



where the generalized partial zeta function 
m n gl J (m n zi, m n Zj) 

d ( ab ) ( \\ ' 

Q u Qv UZ b J | u — n i n 

is introduced with the notation 

a s = (Di + s - l)/2. 



(36) 



(37) 



(38) 



The computation of the VEV of the surface energy- 
momentum tensor requires the analytic continuation of 
the function 3>,(s) to the value s = — 1 (here and below 
| s= _i is understood in the sense of the analytic continu- 
ation), 



=-n«)fcgC J ^(«)| i= _ 1 . 



(39) 



In order to obtain this analytic continuation we will 
follow the procedure multiply used for the evaluation of 
the Casimir energy (see, for instance, The starting 

point is the representation of the function l|37|) in terms 
of the contour integral 



0/3 (s) 



r(-a s )^- s - 

2mT(-s/2) J c 
gj l \uzi,uzj) 

X (ab), 7' 

gl '(uz a ,uz b ) 



duu(u 2 + \ 2 ) c 



(40) 



where C is a closed counterclockwise contour in the com- 
plex u plane enclosing all zeros m n ■ The location of these 
zeros enables one to deform the contour C into a seg- 
ment of the imaginary axis (—iR, iR) and a semicircle 
of radius R, R — > oo, in the right half-plane. We will 
also assume that the branch points are avoided by 

small semicircles in the right half-plane with small 



radius p. Here we assume that there is no zero mode cor- 
responding to m n — 0. The changes introduced by the 
presence of the zero mode will be discussed below in this 
section. For negative Re s with sufficiently large absolute 
value the integral over the large semicircle in Eq. (|40|) 
tends to zero in the limit R — > oo , and the expression on 
the right can be transformed to 



0/3 (s) 



E 

.a=± 



du 



sin 7ra, 



r(-q«) 
r (-j)M*- 1 

gj'\uzi,uzj) 

(ab) / \ 

gl '(uz a ,uz b ) 

v /„2 \2\a 3 Gy\uZi,UZj) 
Gl '(UZ a ,UZb) 



{u 2 + \lY 
2ni 

du u 



TV 



A/3+p 



(41) 



For —D\ — 1 < Res < —T>\ the second integral on the 
right is finite in the limit p — ► 0, whereas the first integral 
vanishes as p as+1 . As a result, by using the formula 
r(— a s ) sin7ra s = — ir/T(a s + 1), in this strip of complex 
plane s we have the following integral representation 



cm = C\s) 



where 



v -f)r(a s + i) 

duu{u 2 - \ 2 p) as Q, jv {uz a ,uz b ), (42) 



n^p- 8 - 1 



T{-s/2)T{a s + 1) 



F^\uzj) ' 



(43) 



In the last formula we use the notation F = K for j = a 
and F = I for j = fe. The contribution of the second 
term on the right of Eq. I|42l) is finite at s = — 1 and 
vanishes in the limits z a — > or z — > oo. The first 
term corresponds to the contribution of a single brane 
at z = Zj when the second brane is absent. The surface 
energy density corresponding to this term is located on 
the surface y = a + for the brane at y = a and on the 
surface y = b— for the brane at y = b. To distinguish on 
which side of the brane is located the corresponding term 
we use the superscript J=L for the left side and J=R 
for the right side. As we consider the region between 
the branes, in formula l|43|l J=R for j = a and J=L for 
j = b. The further analytic continuation is needed for the 
function ^j^(s) only and this is done in the next section. 

In the discussion above we have assumed that there 
is no zero mode corresponding to m n = 0. For spe- 
cial values of the parameters of the model, satisfying the 
condition 



(Aj + B a v)(A b — B b v) ( — 

z b 



(A a -B a is){A b + B b v), 
(44) 
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the value m n = is a solution to the eigenvalue equation 
(|23[) and the zero mode is present. We will denote this 
mode by mo = 0. The corresponding contribution to the 
VEV of the surface energy-momentum tensor and to the 
zeta function defined above is taken into account if wc 
assume that the summation over n in the corresponding 
formulas includes the summand with n = as well. The 
expression for this term is obtained from the formula for 
general n taking the limit m n — > 0. In particular, the 
contribution of the zero mode to the partial zeta function 
(|3*7jl is given by 



T(-a s ) 



\2a 
A /3 



•gl l) (0,0) 



r(- s /2)M s+1 # 



ifjxg v '[uz a ,uz b )\ u= o 



(45) 



As before, the contribution of the nonzero modes, can 
be presented as the contour integral given by Eq. (|40|l . 
but now the contour C contains an additional semicircle 
Cp°^ of small radius which avoids the origin in the right 
half-plane. Consequently, this contribution is given by 
the right hand side of Eq. I|41|) plus the integral over 

the contour C^p \ Evaluating the latter integral we can 
see that it cancels the contribution coming from the zero 
mode presented by Eq. Hence, we conclude that the 

integral representation for the partial zeta function given 
by Eqs. (|4*2|) . t4l?|> is valid for the case of the presence of 
the zero mode as well. 



III. ENERGY DENSITY ON A SINGLE BRANE 

In this section we will consider the geometry of a sin- 
gle brane placed at y = j. The orbifolded version of 
this model corresponds to the higher dimensional gen- 
eralization of the Randall-Sundrum 1-brane model with 
the brane location at the orbifold fixed point y = 0. The 
corresponding partial zeta function is given by Eq. (|43|l . 
Now in this formula J — L, R for the left and right sides of 
the brane, respectively, and F = K for J = R and F = I 
for J = L. In addition, the replacement rft' — > n^' 
should be done with rfi^ = 1 and rfi^ = —1. The in- 
tegral representation (|43[) for a single brane partial zeta 



function is valid in the strip — D± — 1 < Res < — D\ 
and under the assumption that the function F„ (u) has 
no real zeros. For the analytic continuation to s = — 1 
we employ the asymptotic expansions of the modified 
Bessel functions for large values of the argument |57| . 
For Bj from these expansions one has 



1 



i=i 



(46) 



where the coefficients ul .'^{y) are combinations of the 
corresponding coefficients in the expansions for the func- 
tions F v (u) and F' v (v). Note that one has the relation 



.XK,j) 



vi -=(-1)S (/j) , (47) 
assuming that the coefficients in the boundary conditions 
are the same for both sides of the brane. For the nonzero 
modes along the internal space S we subtract and add 
to the integrand in l|43|l the N leading terms of the cor- 
responding asymptotic expansion and exactly integrate 
the asymptotic part: 



r(-s/2)/i s+1 

( F u {u) 



.DC 



du 



(u 2 



\2 „2\u s 

-V.7 ) 



,( J 

a 



N (F,j) 1 



^ 2Bi 

i=i 3 



vDi+s-l + l 



ra/2) 



a, - 1 



(48) 



For the zero mode we first separate the integral over the 
interval (0, 1) and apply the described procedure to the 

integral over (I, oo). By using these formulas for rf?(s), 
the corresponding function 

= Sl^(*)l a <i?00. ( 49 ) 



is written in the form 



$f (s) = 



(4 7 r)( 1 -^)/ 2 n( J )^f 2 f _ 



r(-|)r( 



duu Dl+s 



+ / duu(u^ - ^$Zj) c 

J U/3 



B 



F„{u) 



N F i N 

y v i 



JV 

E 



F v (u) 
F { J\u) 

Dx+s-l+l 



2 r(a ' + 1) g r ( z/2) T {-2 



a s -l]Cs( --a s -l,X 



(50) 



where up — A^Zj+^oa^ and a s is defined by formula (|38|l . function associated with the massless laplacian defined 
In Eq. (|5n)) we have introduced the local spectral zeta 
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on the internal subspace S: 



(51) 



where the prime on the summation sign means that the 
zero mode should be omitted. Both integrals in Eq. H5U|) 
are finite at s — —1 for N ^ D\ — 1. For large values 
A/3 the second integral behaves as \p 1+s ~ N and the series 
over P in Eq. (|50|l is convergent at s = — 1 for N > D — l. 
For these values N the poles at s = — 1 are contained only 
in the last two terms on the right. 

The zero mode part has a simple pole at s = — 1 pre- 
sented by the summand I = D — 1 of the second sum in 
figure braces. The pole part corresponding to the nonzero 
modes is extracted from the pole structure of the local 
zeta function (|51|l . The latter is given by the formula 



r(z)c s (z,x)| 2 



Cd 2 /2- p (X) 



z — p 



n p (x) 



(52) 



where p is a half integer, the coefficients Cd 2 /2- p (X) 
are related to the Seeley-DeWitt or heat kernel coeffi- 
cients for the corresponding non-minimal laplacian, and 
the dots denote the terms vanishing at z = p. In the 
way similar to that used in Ref . [13] , it can be seen that 
the coefficients C P (X) are related to the corresponding 
coefficients C p (X,m) for the massive zeta function 



Cs(s,X;m) 



v p 



m 2 ) 8 



by the formula 



C P (X) = C p (X,0) - \MX)\ 2 S p ,d 2 /2- (54) 



The VEV of the energy density on a single brane is de- 
rived from 



£ f = -»WfcgG J *?>(«)| J= _ 1 . 



(55) 



By using relation (|52[1 . the energy density is written as a 
sum of pole and finite parts: 



(J) _ (J) (J) 



(56) 



Laurent-expanding the expression on the right of Eq. I|5()|l 
near s = — 1, one finds 



.(J) 



oh D r JLJ F 'iK.D-i 



(^^^(s + i) ^ r(i/2) 

x [C (D _ l)/2 (X) + \MX)\ 2 S lDl ] (57) 



(53) for the pole part, and 



4? = 2k°z^C J P Dl J2\MX)\ 2 



SoxpBj I duu 



D-l F A U ) 



+ r duu{u 2 - xy^ 2 - 1 [B j 




N 



-kgzP* c 3 p Dl \Mx)\ 2 { 



i 2v 



1=1 



D 1 I 



h D C N (F,j) 



(47r) B i/ 2 f-^ 



1=1 



21n/i-V( 2 



( j y) - t> (J 



+ fy;-r>i)/2p0 



r 



(58) 



for the finite part, with VK^) being the diagamma func- 
tion. In this formula the prime on the summation sign 
means that the term with I = D\ should be omitted and 
it is understood that C P (X) = for p < 0. In the pole 
part the second term in the square brackets comes from 
the zero mode along £ and this term is cancelled by the 
delta term on the right of Eq. I|54|l . For a one parame- 
ter internal manifold £ with the length scale L one has 
A/5 ~ 1/L and ipp(X) ~ l/L Da . In this case instead of the 



zeta function (|51fl we could introduce the dimensionless 
function Cs(s, X) = L D2 ~ 2s (z(s,X). The corresponding 
dimensionless coefficients Cd 2 /2- p {X) and £l p (X) in the 
formula analog to Eq. (|52|l , are related to the coefficients 
in Eq. I|52|> by the formulas 



Cd 2 /2- p (X) - L 2p D2 C D2 / 2 -p(X), 

fl p (X) = L 2 p- d > \n p {X) + 2C D2/2 ^ p {X) InL 



(59) 
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and do not depend on L. Now we see that the term with 
In L is combined with the term In \x of Eq. (|58|l in the 
form ln(/xi). 

The renormalization of the surface energy density can 
be done modifying the procedure used previously for the 
renormalization of the Casimir energy in the Randall- 
Sundrum model 0, IT^ . IT^ | and in its higher-dimensional 
generalizations with compact internal spaces |26l |27| . 
The form of the counterterms needed for the renormal- 
ization is determined by the pole part of the surface en- 
ergy density given by Eq. I|57l) . For an internal man- 
ifold with no boundaries, this part has the structure 

J2i=o ^ a< t\ z j I L) 21 ■ By taking into account that the 
intrinsic scalar curvature Rj for the brane at y = j con- 
tains the factor (zj/L) 2 , we see that the pole part can be 
absorbed by adding to the brane action counterterms of 
the form 



[(D-l)/2] 

hi £ b^R^ 

1=0 



(60) 



where the square brackets in the upper limit of summa- 
tion mean the integer part of the enclosed expression. 
By taking into account that there is the freedom to per- 
form finite renormalizations, we see that the renormal- 
ized surface energy density on a single brane is given by 
the formula 



(J,ren) _ (J) 
£ ~ £ 3,i 



[(-D-l)/2] 

E 



(61) 



The coefficients c> 8 ' in the finite renormalization terms 
are not computable within the framework of the model 
under consideration and their values should be fixed by 
additional renormalization conditions. 

The total surface energy density for a single brane at 
y = j is obtained by summing the contributions from the 
left and right sides: 



_(LR) 



+ e 



(R) 



(62) 



In formulas (|5T|) , JSSJ we should take F = I for J = L 
and F = K for J = R. Now we see that, assuming the 
same boundary conditions on both sides of the brane, the 
coefficients C P (X, 0) enter into the sum of pole terms in 
the form 



E 

i=i 



r(i) 



21 



c 



D/2-l 



(X,0). 



(63) 



If the internal manifold contains no boundaries and D is 
an odd number, one has Cd/ 2 -i{X, 0) = and, hence, 
the pole parts coming from the left and right sides cancel 
out. In this case for the total surface energy density one 
obtains the formula 



_(LR) 



/oo 
duu(u 2 -X 2 z 2 ) Dl/2 ^\ 
-.a 



du u 



D-l 



/„(«) , K v {u) 



li j \u) Kb J '(u) 



0) / 



I v (u) , K v {u) 



[n/2] {I j) • 



D?\u) Ri j \u) 



,21 



1 = 1 



-2\MX)\ 2 E 



U 2l 



1=1 



Di - 21 



(64) 



r 



where the coefficients v\ \ T are defined by relation (|46|l . 
Note that this quantity does not depend on the renor- 
malization scale (i. In the orbifolded version of the 
model with a single brane at y = 0, which corresponds 
to the higher dimensional generalization of the Randall- 
Sundrum 1-brane model, the bulk is symmetric under the 
reflection y — > —y. In this model the surface densities on 
the left and right sides of the brane are the same and co- 
incide with £j R '- In particular, here the abovementioned 
cancellation of the pole parts from left and right sides 
does not take place. 



For a one parameter internal space of size L the surface 
energy density on the brane at y = j is a function on the 
ratio L/zj only. Note that in the case of the AdS bulk 
the corresponding quantity does not depend on the brane 
position. To discuss the physics from the point of view 
of an observer residing on the brane, it is convenient to 
introduce rescaled coordinates 

x' M = e- kDj x M , M = 0, 1, ...,£>- 1. (65) 

With this coordinates the warp factor in the metric is 
equal to 1 on the brane and they are physical coordi- 
nates for an observer on the brane. For this observer the 
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physical size of the subspace E is Lj — Le and the 
corresponding KK masses are rescaled by the warp fac- 
tor: \p = \pe . Now we see that the surface energy 
density is a function on the ratio Lj/(l/kr>) of the phys- 
ical size for the internal space (for an observer residing 
on the brane) to the AdS curvature radius. 

As an application of the general results presented 
above, we can consider a simple example with E = S . In 
this case the bulk corresponds to the AdS,o + i spacetime 
with one compactified dimension X. The corresponding 
normalized eigenfunctions and eigenvalues are as follows 



a 1-ni0XjL 



13 = 0,±1,±2, 



(66) 



where L is the length of the compactified dimension. The 
surface energy density induced on the brane is obtained 
from general formulas by the replacements 



J2\MX)\ 2 ^lf^', D- 2 1. ((,7) 

0=0 



where the prime means that the summand /3 = should 
be taken with the weight 1/2. For the local zeta function 
from Eq. 115 1H one has 



0=1 



-2s 



2L 2 



(27T) 



-Ch(2s), (68) 



where Cr( z ) is the Riemann zeta function. Now the 
only poles of the function T(z)Cs(z, X) are the points 
z = 0,1/2. By using the standard formulas for the 
gamma function and the Riemann zeta function (see, for 
instance, |57J), it can be seen that one has 

C 1/2 (X) = -1/L, C (X) = 1/2^, (69) 

for the residues appearing in i|52|) and 



O (X) = l4^,O iW = 7 + 21n(L/47r) 



L 



20? 



(70) 



2L 2p ~ 1 1 
top(X) = ^2f T (p)U2p), P + 0, -, 

for the finite parts, with 7 being the Euler constant. 



IV. TWO-BRANE GEOMETRY AND INDUCED 
COSMOLOGICAL CONSTANT 

As it has been shown in Section [D] the partial zeta 
function related to the surface energy density on the 
brane at y = j is presented in the form l|42|l . where the 
second term on the right is finite at the physical point 
s = — 1. By taking into account formulas l|36|) . (|39fl . for 



two-brane geometry the VEV of the surface energy den- 
sity on the brane at y = j is presented as the sum 



.0) 



-(J) 



As 



00 



(71) 



du r i 



(72) 



The first term on the right is the energy density induced 
on a single brane when the second brane is absent. The 
second term is induced by the presence of the second 
brane and is given by the formula 

PC 

Aef = 2C 3 n^{k D z ] ) D Bp Dl Y J \MX)\ 2 

Jx t 

x(u 2 - \l) D ^ 2 - l Sl jv {uz a ,uz b ). 

By taking into account relation © between the VEV of 
the field square and the surface energy-momentum ten- 
sor, this formula can aslo be obtained from the last term 
of Eq. H28f) evaluated at the brane z — Zj. As we con- 
sider the region a ^ y ^ b, the energy desnity (|72() is 
located on the surface y = a + for the left brane and on 
the surface y = b — for the right brane. Consequently, 
in formula l|71|l we take J = R for j = a and J = L for 
j = b. The energy densities on the surfaces y = a — and 
y = b + are the same as for the corresponding single 
brane geometry The expression on the right of Eq. H72|) 
is finite for all nonzero distances between the branes and 
is not touched by the renormalization procedure. For a 
given value of the AdS energy scale ku and one param- 
eter manifold E with the length scale L, it is a function 
on the ratios Zb/z a and L/z a . The first ratio is related 
to the proper distance between the branes, 



Zb/z a 



(73) 



and the second one is the ratio of the size of the internal 
space, measured by an observer residing on the brane at 
y = a, to the AdS curvature radius k^ 1 . The expression 
(|72|l for the surface energy density can be presented in 
another equivalent form: 



Ae 



to _ 







du u(u — Aa) 



2 \Di/2-l 



2k°z? +1 Bp Dl C J 
B?(u 2 Zj + v 2 ) — A 2 j 

_d_ 

dz. 



In 



ll a) (uz a )Kl u >(uz b ) 



ii b \uz b )Rr(uz a ) 



(a), 



(74) 



which will be used below in the discussion of the relations 
between the bulk and surface energy densities. 

For the comparison with the case of the bulk spacetime 
AdSd+i when the internal space is absent, it is useful in 
addition to the VEV l|72() to consider the corresponding 
quantity integrated over the subspace E: 



Ae 



(s) 
Dij 



( s ) e -D 2 k D j 







(75) 
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where Ae^g is defined by the relation 



(76) 



Comparing this integrated VEV with the corresponding 
formula from Ref. [351 ]. we see that the contribution of 
the zero KK mode (Xp = 0) in Eq. H75JI differs from the 
VEV of the energy density in the bulk AdSz^+i by the 
order of the modified Bessel functions: for the latter case 
v — > v\ with v\ defined by Eq. I|20[) with the replacement 
D — > D\. Note that for £ ^ Cd+-Di+i one has v ^ v\. In 
particular, this is the case for minimally and conformally 
coupled scalar fields. 

Now we turn to the investigation of the part (|72l) in the 
surface energy density in asymptotic regions of the pa- 
rameters. For large values of AdS radius compared with 
the interbrane distance, ko(b ~ a) <C 1, the main con- 
tribution to the integral on the right of Eq. I|72l) comes 
from large values of uz a ~ [ko(b — a)]^ 1 . Assuming that 
B a /(b — a) and m(b — a) are fixed, we see that the or- 
der of the Bessel modified functions is large. Replacing 
these functions by their uniform asymptotic expansions 
for large values of the order [57]], one obtains 



/c 

(u 2 - ml)^-\M - mu 2 )- 1 



[c a (u)c b (u)e 2ui ~ b - a ) - l] 



du u 



(77) 



where mp 
tation 



i 2 + \ 2 p and we have introduced the no- 



j=a,b. (78) 



U) Bj 
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The expression on the right of Eq. I|77|l is the correspond- 
ing surface energy on the brane in the bulk geometry 

For large KK masses along E, z a Xp ^> 1, Xp ^> 1, we 
can replace the modified Bessel functions by the corre- 
sponding asymptotic expansions for large values of the 
argument. For the contribution of a given KK mode to 
the leading order this gives 

/>oo 

Aeg « 4nWk%zP +1 B]C j l3 Dl J duu 2 



{A 2 ~ Byz 2 )' 1 ^ 2 - \ 2 )^- 

c a (uz a )c b (uz b )e 2u( - z >>- 



Za) 



; (79) 



where 



cj(u)= i ■ mp J , j = a,b. (80) 



A, 



If in addition one has the condition Xp(z b — z a ) 3> 1, the 
dominant contribution into the it-integral comes from the 



lower limit and we have the formula 



Ae 



(s) 



2n^B 2 C ] k°zf +1 
" ^ A) - (XffZjBj) 2 (4tt)^/2 

A ^/2+l e -2A,( 26 - Za ) 



Ca{^0Z a )c b {XfSZ b )(z b - Z a ) D ^/ 2 ' 



(81) 



In particular, for sufficiently small length scale of the 
internal space this formula is valid for all nonzero KK 
masses and the main contribution to the surface densities 
comes from the zero KK mode. In the opposite limit of 
large internal space, to the leading order we obtain the 
corresponding result for parallel branes in AdSo+i bulk 

M 

For large values of the mass with m 3> m(z b — 
z a ) 3> 1, and m 3> Xp, one has v ~ m/ko S> 1. Intro- 
ducing a new integration variable u — vy and using the 
uniform asymptotic expansions for the modified Bessel 
functions, it can be seen that the corresponding asymp- 
totic formula for Aefl is obtained from l|81|l by the re- 



placement Xp — > m. If both masses m and A/3 are of 

(s) 

the same order, the asymptotic behavior for Ae^g is de- 
scribed by Eq. l{8"Tll replacing Xp — > mp. As we could 
expect in this case the induced surface densities are ex- 
ponentially suppressed. 

For small interbrane distances, ko(b~a) <C 1, which is 
equivalent to z b /z a - 1 < 1, the main contribution into 
the integral in Eq. (|72|l comes from large values u and to 
the leading order we obtain formula (|79|l . If in addition 

one has Xp{z b - z„) < 1 or equivalcntly X^\b — a) <C 1, 
we can put in this formula Xp = 0. Assuming (6 — a) <C 
\Bj/Aj\, for Bj ^ to the leading order one finds 



Ae 



c 3 r(Y)C^i-i) 



(s) 



U). 



(82) 

where exj = 1 for |£? a /A a |, |_B(,/^4h| ^> ko(b — a), and 
(Tj = 2 2 -° - 1 for {Bj/Aj] > k D (b- a) and B l /A l = 0, 
with I = b for j = a and I = a for j = b. We see that for 
small interbrane distances the sign of the induced surface 
energy density is determined by the coefficient Cj and 
this sign is different for two cases of Uj. 

Now we consider the limit Xpz b 3> 1 assuming that 
Xpz a < 1. Using the asymptotic formulas for the Bessel 
modified functions containing in the argument z b , for the 
contribution of a given nonzero KK mode we find the 
following results 



Ae 



.(b) 
-a/3 

(b) 
hti 



k°z°B 2 a C a (Xp/z b ) D ^ 2 e- 2 ^ 
^ D ^- l c b (Xpz b )K^\Xpz a ) 

kZzF>^B 2 C b A^ 2+ V 2 ^ 



(83) 



2 D 1 -l 7r Z5 1 /2-l z ^i/2 (A h + XpZ b B b y 

li a) {Xpz a ) 



Ki a) {Xpz a ) 



(84) 
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This limit corresponds to the interbrane distances much 
larger compared with the AdS curvature radius and with 
the inverse KK masses measured by an observer on the 
left brane: (6— a) » 1/fco, 1/Ajg . For a single parameter 
manifold £ with length scale L and (b — a) L a these 
conditions are satisfied for all nonzero KK modes. 

In the limit z a X/3 <C 1 for fixed Zb\p, by using the 
asymptotic formulas for the modified Bessel functions for 
small values of the argument and assuming \A a \ ^ \B a \v, 
one finds 



Ae (s) 



2 2 ^T 2 {v){A a -vB a ) 

v /,,2 x2\Di/2-l ^ \uz b ) 
Ih '(UZ b ) 



du u 



2v+l 



Ae 



c a {y) \Z b 



(85) 
(86) 



where we have introduced the notation 



.(6) = 4BjC b f3 Dl 
JvP 2 2 »vT 2 (y) 



du 



0Zb 



ll b)2 (u) 



2 2\Dx/2-\ 



(87) 

For \A a \ = \B a \v we should take into account the next 
terms in the corresponding expansions of the modified 
Bessel functions. The integral in Eq. I|85|) is negative 
for small values of the ratio A b /B b and is positive for 
large values of this ratio. As it follows from Eq. (|86|l . 
for large interbtane separations the sign of the quantity 
Ae^ is determined by the combination (B 2 v 2 — A 2 )Cb of 
the coefficients in the boundary conditions. In the limit 
under consideration the KK masses measured by an ob- 
server on the brane at y = a are much less than the AdS 
energy scale, Al -C &d, and the interbrane distance is 
much larger than the AdS curvature radius. In particu- 
lar, substituting A/3 = 0, from these formulas we obtain 
the asymptotic behavior for the contribution of the zero 
mode to the surface energy density induced by the sec- 
ond brane in the limit z a jz b <C 1. Now combining the 
corresponding result with formulas i)83|) , H84fl , we see that 
under the conditions (b — a) l/ku,L a , and L a kjj > 1, 
the contribution of the nonzero KK modes along £ is 
suppressed with respect to the contribution of the zero 

mode by the factor (z b / z a ) 2l/+Dl/2+ ^ exp(-2A /3 z fc ). 

From the analysis given above it follows that in the 
limit when the right brane tends to the AdS hori- 
zon, Zb — > 00, the energy density Ae^i vanishes as 
e~ 2X f Zb I z® 1 ^ 2 for the nonzero KK mode along £ and as 
for the zero mode. The energy density on the 



vanishes as z 



D2 + D1/2+I 2A 



e 

2;/ 



\pz b f or the 



right brane, Ae^ 

nonzero KK mode and behaves like z b U2 ^ 2u for the zero 
mode. In the limit when the left brane tends to the AdS 
boundary, z a — > 0, the contribution of a given KK mode 

vanishes as z^ +2u for Ae^ and as z 2v for AeS. For 
small values of the AdS curvature radius corresponding 



to strong gravitational fields, assuming \pz a 1 and 
Xp{zb — Za) ^> 1, we can estimate the contribution of 
the nonzero KK modes to the induced energy densities 
by formula 181|) . In particular, for the case of a single 
parameter internal space with the length scale L, under 
the assumed conditions the length scale of the internal 
space measured by an observer on the brane at y = a 
is much smaller compared to the AdS curvature radius, 
L a -C k^ 1 . If L a > kp 1 one has \pz a < 1 and to es- 
timate the contribution of the induced surface densities 
we can use formulas and (|81Jl . and the suppression 
is stronger compared with the previous case. For the 
zero KK mode, under the condition kn(b — a) S> 1 we 
have z a j Zb <C 1 and to the leading order the correspond- 
ing energy densities are described by relations l|85|) and 
(|86|l . From these formulae it follows that the induced en- 
ergy densities integrated over the internal space behave 
as k® 1+1 exp[(DiSj +2i')ko(a — b)] for the brane at y = j 
and are exponentially suppressed. Note that in the model 
without the internal space we have similar behavior with 
v replaced by v\ and for a scalar field with ( < Cd+l>i+i 
the suppression is relatively weaker. 

In Fig. ^and Fig. [21 we have plotted the energy den- 
sities induced on the left and right branes, respectively, 
by the presence of the second brane as functions on the 
size of the internal space £ ~ S 1 and interbrane distance 
for a D — 5 minimally coupled massless scalar field with 
the Robin coefficients A a /B a = —1, A b /Bb = 5. The 
corresponding expressions are obtained from the general 
formulas of this section by the subtitutions l(57)l. 



0.05 




0. 8 



FIG. 1: Surface energy density, Aei s /kjj, induced on the 
brane at y = a as a function of z a /zi, and L/zt for a, D — 5 
minimally coupled massless scalar field in the model with E = 
S and with the Robin coefficients A a /B a = — 1, At/Bb = 5. 



Introducing the rescaled coordinates defined by Eq. 
()65p. after the Kaluza-Klein reduction of the higher di- 
mensional Hilbert action, by the way similar to that in 
the Randall-Sundrum braneworld, it can be seen that 
effective Z?i-dimensional Newton's constant Goxj mea- 
sured by an observer on the brane at y — j and the fun- 
damental (D + l)-dimensional Newton's constant Gd+i 
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0.3 



FIG. 2: The same as in Fig. 0for the right brane. 



are related by the formula 



G Dl] - 



(D - 2)k D G D+1 
\ e ( D - 2 ) k °( b ~ a ) - 



1] 



= (D-2)fe D (b-j) 



where 



-D 2 k D j J 



(88) 



(89) 



is the volume of the internal space measured by the same 
observer. In the orbifolded version of the model an addi- 
tional factor 2 appears in the denominator of the expres- 
sion on the right. Formula 188|) explicitly shows two pos- 
sibilities for the hierarchy generation by the redshift and 
large volume effects. Note that the ratio of the Newton 
constants on the branes, G Dlb /G Dia = e ( 2 -Di)k D (b-a) ^ ig 
the same as in the model without an internal space. For 
large interbrane distances one has Gn ia ~ fccGc+i/Vso, 



G 



Dib 



knG 



D<-*D+ie 



(2-D)k D (b-a) 



/Vzb, and the gravita- 
tional interactions on the brane y — b are exponen- 
tially suppressed. This feature is used in the Randall- 
Sundrum model to address the hierarchy problem. As 
we will see below this mechanism also allows to obtain 
a naturally small cosmological constant generated by the 
vacuum quantum fluctuations (for the discussion of the 
cosmological constant problem within the framework of 
braneworld models see references given in j35j). 

As we have already mentioned, surface energy density 
(|72|l corresponds to the gravitational source of the cos- 
mological constant type induced on the brane at y = j 
by the presence of the second brane. For an observer 
living on the brane at y = j the corresponding effective 
Di-dimensional cosmological constant is determined by 
the relation 



Ap x j = 8ttG 



Dij^Dij 



M 



(90) 



Dij 



where Mr> x j is the Di-dimensional effective Planck mass 

(s) 

scale for the same observer and Ae^, ■ is defined by 



Eq. J7SJ). Denoting by M D+1 the fundamental (D + 1)- 
dimensional Planck mass, G D+1 = M%£, from Eq. ggj 
one has the following relation 



M 



Dij 



M 



D+l 



(zb/z a ) 



D-2 



-1 Vsj 



MSlW (91) 



(D-2)(z b /z j ) D ~* k 



D 



for the ratio of the effective and fundamental Planck 
scales. From the asymptotic analysis for the induced 
vacuum densities given above it follows that for large 
interbrane distances, z b /z a 3> 1, the contribution of the 
modes with \pz b ^ 1, Xpz a < 1 is suppressed by the fac- 
tor (z b /z a ) 2u+Dl/2+ ^ exp(-2A /3 z 6 ) with respect to the 
contribution of the zero mode. For the contribution of 
the modes with z a Xp <C 1, z b \p < 1 from Eqs. JSSJl, 

dHl one has Aeg - k%zj*(z a /z b ) 2u+Di '~ D '- By us- 
ing these relations, one obtains the following estimate 
for the ratio of the induced cosmological constant l|9U|) 
to the corresponding Planck scale quantity in the brane 
universe: 



3 ZkGd^M** 



2v+D l 8° 



M, 



Dij 



Di 



(92) 



By taking into account relation (|91|) . this can also be 
written in the form 



k Di-l 

V*jM%-l 



Di/(Di-2) 



x exp 



k D {a-b)[2v + D 1 + 



D2Dl rb 

Di-2 j 



.(93) 



For the model without an internal space this ratio is of 
the same order of magnitude for both branes. 

In the higher dimensional version of the Randall- 
Sundrum braneworld the brane at z = z b corresponds 
to the visible brane. For large interbrane distances, by 
taking into account Eq. (|86|) , for the ratio of the induced 
cosmological constant (|9UII to the Planck scale quantity 
in the corresponding brane universe one obtains 



1 



CD 



c a {v) \M Dlb 



Di 



Zb 



2// 







(94) 



where the function fj^ is defined by Eq. (f8Tj) . In Fig. 
we have plotted the coefficient in this formula as a 
function of B b and L/z b for aD = 5 minimally coupled 
massless scalar field in the model with the internal space 
S = S 1 and A b — 1. We recall that L/z b = koL b , where 
L b is the physical size of the internal space for an observer 
living on the brane y = b. 

Using relation Ij91|l with j = b, we can express the 
corresponding interbrane distance in terms of the ratio 
of the Planck scales 



Zb 



Kit 



1 Vsb ( Mn+i 
k D \ M Dlb 



1 l/P-2) 



(95) 
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FIG. 3: Coefficient f^J in formula 194H as a function on 
Bb and L/zb for a D — 5 minimally coupled massless scalar 
field in the model with £ = S 1 and At = 1. 



Substituting this into Eq. (|94() . for the ratio of the cos- 
mological constant on the brane at j = b to the corre- 
sponding Planck scale quantity one finds 

with i> = 2v/(D — 2). The higher dimensional Planck 
mass Mfl + i and AdS inverse radius fco are two funda- 
mental energy scales in the theory which in the Randall- 
Sundrum model are usually assumed to be of the same 
order, ko ~ Mp+i (see, e.g., [lj). In this case one ob- 
tains the induced cosmological constant which is expo- 
nentially suppressed compared with the corresponding 
Planck scale quantity on the visible brane. In the model 
with D x = 4, k D ~ Md+i ~ 1 TeV, M Dlb = M P1 ~ 10 16 
TeV, assuming that the compactification scale on the 
visible brane is close to the fundamental Planck scale, 
V-£bM®^ 1 ~ 1, for the ratio of the induced cosmolog- 
ical constant to the Planck scale quantity on the visi- 
ble brane we find the estimate h b ~ io~ 32 ( 2+I/ ). From 
(|95|l one has /c_d(6 — a ) ~ 74/(£>2 + 2) and the corre- 
sponding interbranes distances generating the required 
hierarchy between the electroweak and Planck scales are 
smaller than those for the model without an internal 
space. In the model proposed in Ref. [26j. a separa- 
tion between the fundamental Planck scale and curva- 
ture scale is assumed: ko ~ Mn+iz a / z b ~ 1 TeV. Under 
the assumption V^bMj^+i ~ •"■) m this model we have 
fc 6 - i -64[i+>V(c+i)] and _ a ) ~ 74/{D 2 + 3). 



V. TOTAL VACUUM ENERGY AND THE 
ENERGY BALANCE 

On background of manifolds with boundaries the total 
vacuum energy is splitted into bulk and boundary parts. 
In the region between two branes the bulk energy per 
unit coordinate volume in the D\ -dimensional subspace 
is obtained by the integration of the [j-component of the 
volume energy-momentum tensor over this region: 

E (v) = J d° 2 Xdy v^j(0|T o (v)O |0). (97) 

The surface energy per unit coordinate volume in the 
Di-dimensional subspace, E^ s \ is related to the surface 
densities evaluated in previous sections by the formula 

** = li' (98) 

j=a,b J 

Now by making use of the formula for the volume energy- 
momentum tensor from Ref. 01 an d the formula l|34|) 
for the surface densities, it can be seen that the following 
formal relation takes place for the unrenormalized VEVs: 

E = E^+E ( - S \ (99) 

where 

s = ^/r£pr5:£ (fc2 + m " + A ' )1/2 ' (100) 

J ^ ' p n=l 

is the total vacuum energy per unit coordinate volume 
of the I?! -dimensional subspace, evaluated as the sum of 
zero-point energies of elementary oscillators. 

The total vacuum energy within the framework of the 
Randall- Sundrum braneworld is evaluated in Refs. 0, 
IT^ | by the dimensional regularization method and in 
Ref. |12| by the zeta function technique. Refs. 0,0,0 
consider the case of a minimally coupled scalar field in 
D = 4, and the case of arbitrary ( and D with zero 
mass terms c a and c& is calculated in Ref. 0|. For the 
orbifolded version of the model under consideration with 
Di = 4 and zero mass terms on the branes, the vacuum 
energy is investigated in j2(J by using the dimensional 
regularization. Here we briefly outline the zeta function 
approach in the general case. 

We consider the zeta function related to the vacuum 
energy l|100|) : 

CM = m s+1 / tJ^t E X> 2 + < + ^r" a > 

where, as before, the parameter \x with dimension of mass 
is introduced by dimensional reasons. The vacuum en- 
ergy in the region between the branes is obtained by the 
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analytic continuation of this zeta function to the value 
s = -1: 



E = -C(*)| s =-i. 



(102) 



After the evaluation of the integral over k one obtains 
the formula 



C(«) 



s+1 



r( Sl /2) 



(4^-)(Ci-i)/2 r( s / 2 ) 



(103) 



with the partial zeta function 

oo 

(p( Sl ) = EK + ^P l/2 , si = a + 1 - £>i. (104) 

71=1 

To evaluate the vacuum energy we need to perform the 
analytic continuation of the function l|103|) to the neigh- 
borhood of s = — 1. This corresponds to the analytic 
continuation for (,ft{ s i) to the point s\ = —D\. For this 
we present the partial zeta function i|104|) as the contour 
integral 



C/3(«)=2^ / du(u 2 + \l)- s / 2 ^\ngi ab \uz a ,uz b ), 

(105) 

where the integration contour is the same as in formula 
(|40p. Now by making use of the standard properties of 
the Bessel functions, we see that the parts of the integrals 
over (0, cancel and we find the following integral 

representation 



C (s) = isin^ rdutf-\%)-/* 

* 2 J Ah 



x — lnG^(uz a ,uz b ). 
du 



(106) 



For the further discussion it is convenient to write the 
partial zeta function in the decomposed form 



j=a,b 



(107) 



with the separate terms 
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* 2 A fl 
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d i 

x — In 

du 



ll a \uz a )Kl°>(uz b ) 



1 - 
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Kl a \uz a )rt >(uz b ) 
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(108) 



,(109) 



where F — K for j 
A(p (s) is finite at s 



a and F = I for j = b. The part 
— D\ and the analytic continuation 



is necessary for the parts (s), j — a, b. The latter are 
the partial zeta functions for the geometry of a single 
brane in the regions y a and y ^ b, respectively. On 
the base of Eq. (|107|l similar decomposition can be given 
for the total zeta function 



cw= E c (i) (s)+Ac( S ). 



(110) 



where separate parts are related to the corresponding 
partial zeta functions by the relations similar to Eq. 
(I103|l with si defined by Eq. I|104|) . In deriving the inte- 
gral representation (|1U6[) we have assumed that there is 
no zero mode corresponding to m n — 0. The contribu- 
tion of the possible zero mode to the partial zeta function 
Cff(s) is \p S ■ As regards to the nonzero modes, their con- 
tribution, as before, is given by the contour integral on 
the right of Eq. I|105|) , but now the point u = has to be 
avoided by the small semicircle in the right half-plane. 
The integrals over the parts of the imaginary axis are 
transformed to the expression given by the right hand 
side of Eq. (|106fl . whereas the integral over the small 
semicircle is equal to ~A^ S and cancels the contribution 
from the zero mode. Hence, the integral representation 
(|106fl for the partial zeta function is valid in the case of 
presence of the zero mode as well. 

In order to obtain an analytic continuation for the 
function C^( s ); m the- integrand for the non-zero modes 
we subtract and add the term which is obtained replac- 
ing the function F^\uzj) by the first N terms of the 
corresponding asymptotic expansion for large values u. 
In the zero mode part we split the integral in Eq. H108|l 
into the integrals over (0,1) and (l,oo). The contribu- 
tion of the first integral to the zeta function is finite at 
s = —D\ and the part with the second integral can be 
treated in the way similar to that used for the non-zero 
modes. After the explicit integration of the asymptotic 
part and introducing the notation 



£(, Fj) (u) = -n 



0) 



2q F e r ' 



u Bi 



-F^(u), 



(111) 



with qj — it and qx = 1/tt, for the corresponding zeta 
function the following formula is obtained 
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/iS +l (47r )(l-U 1 )/2 z « 

r(s/2)r(l-ai/2) 
d 



L 



du 

N 

1=1 



£ in S<™(») -E 



iV (Fj) ' 



W (Fj) ' 

si +Z 



2 ' -x *i 



(112) 



where, as before, u^g = XpZj + Sqx^, and we have used 
the relation sin(7r,z)r(l — z) = tt/T(z) for the gamma 
function. In formula 1)112(1 . 



20? 



~(FJ) 



(F,i) 



r(//2) 



, 1 = 1,2,..., 



(113) 



with 



,(Fj) 



being the coefficients in the asymptotic ex- 



pansion of the function ln[Ey i? '^('u)] for large values of 
the argument: 

00 (Fj) 



and we have defined the global zeta function 



-2s 
/3 ■ 



(115) 



employ the global analog of formula l|52[) : 



r(«)CE(«)l* 



c 



D2/2-P 

z — p 



i2 m 



(116) 



As it has been shown in Ref. [27|, the coefficients C p 
in this formula are related to the integrated Seeley- 
DeWitt coefficients C p (m) for the massive zeta function 



Cs(s;m) = EaOfl 



by the formula 



C p — C p (0) — 5 p ,d 2 / 2 , 



(117) 



which is obtained from (|54|l by the integration over the 
internal space. For the model with the internal space E = 
S 1 the zeta function £s (z) is related to the corresponding 
local zeta function discussed in Section lTTTl bv the formula 
(s(z) = LC,y,{z 1 X) and the expressions for the coefficients 
Cd 2 /2- p and ^ p are directly obtained from fflfy . JJJJ. 

By making use of formula (|116|) . near the point s = — 1 
the zeta function Q^' (s) is presented in the form of the 
sum of pole and finite parts: 



Note that one has the relation wi K ^' = (— l) z io! ; , as- 
suming that the coefficients in the boundary conditions 
are the same, and these coefficients are related to the co- 
efficients in the similar expansions for the modified Bessel 
functions (see, for instance, |H3)- For N > D\ both in- 
tegrals on the right of formula l|112|) are finite at the 
physical point s = — 1. For large values of Xp the second 
integral behaves as (XpZj)^ N ~ Sl and, hence, the series 
over f3 is convergent at the point s\ = — D\ ii N > D . 
As a result, for these values of TV the only poles in the 
expression l|112[) are those for the function T(z)(^(z) and 
the pole in the last sum corresponding to the term with 
I = D\. Note that taking N = 00 we obtain the expan- 
sion of the zeta function over l/zj. 

In order to separate the pole and finite parts of the 
zeta function C^( s ) a t the physical point s — —1, we 



Ao) 

c« ) (-)i^- 1 = M+^. 



s+l 



with the residue 



N 



~ (Fj) 



( 47r )D 1 /2 



(118) 



(119) 



In formula l|119|) the part with the first term in the square 
brackets comes from the nonzero KK modes along E and 
the part with the second term comes from the zero mode. 
Now by taking into account formula 1)1171) . we see that 
the zero mode part is cancelled by the delta term in Eq. 
HI 171 ■ For the finite part of the zeta function we obtain 
the formula 
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i=-i J 



2C* (jD _ ;)/2 ( ln/i - ( -- ) ) + Sl(i- Dl )/2 



(120) 



where the prime on the summation sign means that the 
term / = D\ should be omitted. Note that taking in this 
formula N = oo wc obtain the expansion of the finite 
part over 1/zj. For an internal space with the length 
scale L we can introduce the dimensionless zeta function 
Cs(s) = £~ 2s Cs( s ) with the coefficients and & P 

defined by the relation similar to Eq. (|1 . Now it can 
be seen that in terms of these coefficients the last term 
on the right of formula (|120fl is rewritten as 



T-Di N , 
L V w {F ' j) ' 



(An) 
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2C, 
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(121) 



(?) 

Z i 



This explicitly shows that the combination Q 
function on the ratio L/zj only. 

Similar to the case of the surface energy, for a single 
brane at y = j we denote by E^ the total vacuum energy 
in the region y ^ j for J = R and in the region y ^ j 
for J = L. This energy is related to the zeta function 
considered above by the formula 



E 



(J) _ 



(122) 



where in the formulas for the function C ( s ) it should 
be taken F = I for J = L and F = K for J = R. In 
accordance with Eq. i|118|) the regularized vacuum en- 
ergy contains pole and finite parts. The structure of 

( 7*1 

the residue (,_{ determines the form of the counterterms 
which should be added to the action for the renormaliza- 
tion of the vacuum energy. From formula H119(l it follows 



that for an internal space without boundaries the pole 
terms can be absorbed by adding to the brane action the 
counterterm 



\{D+l)/2] 

\h\ aiR\, 



1=0 



(123) 



where Rj is the intrinsic curvature scalar for the brane at 
y = j. Note that this counterterm has the same structure 
as that for the surface energy. The only difference is 
the additional summand corresponding to the term I = 
[(D + l)/2]. Hence, by adding the counterterm l|123|) 
we can absorb the pole terms for both total and surface 
energies. Now by taking into account the possibility for 
the appearance of the finite renormalization terms, the 
renormalized vacuum energy for the geometry of a single 
brane is presented in the form 



£,(J,rcn) _ ^(J) 



E 

1=0 



ci(zj/Lf 



(124) 



where E$ — Co" an d the coefficients q in the finite 
renormalization terms should be fixed by imposing renor- 
malization conditions. 

As in the case of the surface energies, now we see that 
for the internal manifold with no boundaries in the cal- 
culation of the total vacuum energy for a single brane at 
y = j, including the contributions from the left and right 
regions, the pole parts of the energies cancel out for odd 
values D (assuming that the coefficients in the boundary 
conditions JHJl on the right and left surfaces are the same) 
and we obtain a finite result. For this energy one receives 
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Ej = 
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(125) 



with M > (D — l)/2. Note that the energy per unit 
physical volume in the Di-dimensional subspace is de- 
termined as {kjjZj) Dl and is a function on the ratio 
L/zj only. 

On the base of the analysis given above, for the geome- 
try of two branes the total energy in the region a ^ y ^ b 
is presented in the form 



E = E^ + E, 



(L) 



AE, 



(126) 



where the interference term AE in this formula is finite 
for all nonzero values of the interbrane separation and is 
obtained directly from the part A£ (s) substituting s = 
— 1. After the integration by parts this term is presented 
in the form 



AE 
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du u(u — Xa) 



2 \Di/2-l 



x In 
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(uz a )K„ \uzb) 
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(127) 



and is not affected by finite renormalizations. 

Now we briefly discuss the behavior of the interfer- 
ence part in the vacuum energy in the asymptotic regions 
of the parameters. For large values of the AdS radius 
compared with the interbrane distance, ko{b — a) <C 1, 
by making use the uniform asymptotic expansions for 
the modified Bessel functions, to the leading order one 

finds the result for two parallel branes on the bulk 
R ( Dl -i,i) x s . 
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x In 



du u(u 2 
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c a (u)c b (u) 



(128) 



where mp — ^Jm 2 + A 2 and Cj(u) is defined by formula 

(|78|l . For large KK masses along the internal space, 
z a Xp 2> 1, for the contribution of the mode with a given 
(3 one has 



AE B 



duu(u 2 - X}) ^ 2 - 1 



e -2u(z h -z a ) 

x In 1 — 



C a (uZ a )Cb(uZb) 



(129) 



with Cj(u) defined by Eq. (JSDJ. This contribution is 
exponentially suppressed if Xp(zb — z a ) 1: 
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2(47r) z V 2 \zb~z a ) c a (Xi3Z a )cb(Xf}Zb)' 

(130) 

For small interbrane distances, satisfying the conditions 
k]j(b — a) 1 and Xp{zb — z a ) <C 1 to the leading order 
we have the formula 



AE, 







■r i) e - D ^ k ° a (Dt 



(4 7 r)( £, i+ 1 )/2(6-a) £, i 



1 



(131) 



where a = 1 for c a (Xpz a )cb(XpZb) > and a = 2~ Dl — 1 
for c a (X{3Z a )cb(Xf3Zb) < 0. Noting that the renormalizcd 
vacuum energies for single branes are finite in the limit 
a — > 6, we see that for sufficiently small interbrane dis- 
tances the totall vacuum energy is dominated by the in- 
terference part. In particular, it follows from here that 
for the case a = 1 any fixation of the interbrane distance 
can be stable only locally. 

In the limit Xpz b 3> 1 assuming that Xpz a < 1, for the 
contribution of a given KK mode one has the estimate 



AE. 



'0 



V 47re6/ c b (X(3z b ) Ki a) {X p z a y 



(132) 



and the contribution from this part of KK spectrum is 
exponentially small. For Xpz a <C 1 and for fixed Xpz b we 
find 



)l-2v 



AE B 



2v r°° 



■jT 2 (v) Ca {v) Jx 



duu 2v+1 



x(u a -A|)^^j™. (133) 
lb (uz h ) 

From the asymptotic analysis given above it follows that 
when the right brane tends to the AdS horizon, z b — > oo, 
one has the estimate AEp ~ e~ 2X ' iZb / z^ 1 ^ 2 for the 
nonzero KK modes and AEp ~ z~ Dl - 2v for the zero 
KK mode. When the left brane tends to the AdS bound- 
ary, z a — ► 0, the interference part of the vacuum en- 
ergy behaves as AE ~ z 2 ". For small values of the 
AdS curvature radius corresponding to strong gravita- 
tional fields the main contribution comes from the zero 
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mode and the interference part is suppressed by the fac- 
tor e ^k D (a-b) J z D t ^ 

In figures we have presented the dependence of the 
interference part of the vacuum energy on the size of the 
internal space X = S 1 and interbrane distance for afl = 
5 minimally coupled massless scalar field with the Robin 
coefficients A a /B a = -1, A b /B b = 5 (Fig. and B a = 
0, A b /B b = 5 (Fig. [5}. Note that, in accordance with 
the asymptotic analysis, for small interbrane distances 
this part is negative for the first case and positive for the 
second one. In the second example the interference part 
of the vacuum energy has minimum with respect to both 
variables 




FIG. 4: Interference part of the vacuum energy, AE/kp +1 , as 
a function of z a /zb and L/zb for a D — 5 minimally coupled 
massless scalar field in the model with £ = S 1 and with the 
Robin coefficients A a /B a = —1, Ab/Bb = 5. 




0. 7 



0.5 



FIG. 5: The same as in Fig. 21 for the values of the Robin 
coefficients B a — 0, Ab/Bb = 5. 

Now let us check that for the separate parts of the vac- 
uum energy the standard energy balance equation takes 
places. We denote by P the perpendicular vacuum stress 
on the brane integrated over the internal space. This 
stress is determined by the vacuum expectation value of 
the ^-component of the bulk energy-momentum tensor: 



P = - J d D2 X^{Q\T^ )D \0). We expect that in the 
presence of the surface energy the energy balance equa- 
tion will be in the form 

dE = -PdV + E f d S U \ Ef = J d D2 X^ef, 

j—a,b 

(134) 

where V is the (D + l)-volume in the bulk and S (3 > is the 
D-volume on the brane y — j per unit coordinate volume 
in the Di -dimensional subspace, 



V = J dye- DkDV J d D *X^y, 
S U) _ e -Dk D j J d D2 X ^y, j = a, b. 

Combining equations (|134|) - (|136[) . one obtains 
dE n^pti) - Dk D E<f ] 



(135) 
(136) 

(137) 



with Pw) being the perpendicular vacuum stress on the 
brane at y = j. The vacuum stresses normal to the 
branes can be presented in the form p^' = p{ + p 



j = a, b, where p^ is the stress for a single brane at y = j 

(i) 

when the second brane is absent and p)^ t \ is induced by 
the presence of the second brane. The latter determines 
the interaction forces between the branes and is defined 
by the formula given in Ref. ^ij. Here we consider 
the corresponding quantity integrated over the internal 
space: 



(int)> 



p(j) _ 
^(int) 



= ld^x^> t) =k 



/>oo 

(3 Dl / dwi 



x(u 2 - Xl^-^iuz^uzbJf^iuzjUUS) 



with the notation 



u 2 - 



v 2 + 2^B*-D{At-l)A j B j -A 



. (139) 

As in the case of the induced surface energy density, the 
formula for the interaction force can also be presented in 
the form 



P, 



(int) 
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(k D z 3 ) D+1 (3 Dl J2 / duu 

8 J X P 



:(u 2 - X 2 )^' 2 - 1 



2DB 2 C 3 



d i 



li a \uz a )Kr(uz b ) 



B 2 (x 2 z 2 + v*) - A 2 
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(140) 



Now by taking into account expressions l|74|l . (|127|l . 
(114m . it can be explicitly checked that the interference 
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parts in the vacuum energies and effective pressures on 
the branes obey the energy balance equation 

dAE _ nWP$ t) -Dk D AEM 

The second term in the nominator on the right of this 
formula corresponds to the additional pressure acting on 
the brane due to the nonzero external curvature. 



VI. CONCLUSION 

Continuing our previous work |4.ll l44j , in the present 
paper we have investigated the expectation value of the 
surface energy-momentum tensor induced by the vacuum 
fluctuations of a bulk scalar field with an arbitrary curva- 
ture coupling parameter satisfying Robin boundary con- 
ditions on two parallel branes in background spacetime 
AdSuj+i x £ with a warped internal space S. Vacuum 
stresses on the brane are the same for both subspaces and 
the energy-momentum tensor on the brane corresponds 
to the source of the cosmological constant type in the 
brane universe. It is remarkable that the latter prop- 
erty is valid also for the more general model with the 
mertic g^ a instead of rj^ a in line element J2J). As an reg- 
ularization procedure for the surface energy density we 
employ the zeta function technique. By using the residue 
theorem we have constructed an integral representations 
for the partial zeta function corresponding to a given 
KK mode along the internal subspace. This function is 
presented as the sum of single brane and second brane 
induced parts [see formula (|42Jl ]. The latter is finite at 
the physical point and the further analytical continua- 
tion is necessary for the first term only. We have done 
this in Section ITTT1 As the first step we subtract and add 
to the integrand the leading terms of the corresponding 
asymptotic expansion for large values of the argument 
and explicitly integrate the asymptotic part. Further, 
for the regularization of the sum over the modes along 
the internal space we use the local zeta function related 
to these modes. By making use of the formula for the 
pole structure of this function, we have presented the en- 
ergy density on a single brane as the sum of the pole 
and finite parts. These separate parts are determined by 
formulas l|57|) . I|58|l . The pole parts in the surface en- 
ergy density are absorbed by adding to the brane action 
the counterterms having the structure given by Eq. l(B0j l. 
The renormalized energy density on the corresponding 
surface of a single brane is determined by formula (|61ll . 
where the second term on the right presents the finite 
renormalization part. The coefficients in this part can- 
not be determined within the model under consideration 
and their values should be fixed by additional renormal- 
ization conditions which relate them to observables. 

Unlike to the single brane part, the surface energy den- 
sity induced by the presence of the second brane con- 
tains no renormalization ambiguities and is investigated 



in Section IIVI This part is given by formula l|72|l , or 
equivalently by formula (|71J) . We have investigated the 
induced energy density in various asymptotic regions for 
the parameters of the model. In the limit ko —> we 
have obtained the result for the geometry of two bound- 
aries in the bulk ijt^ 1-1 ' 1 ) x S. For the modes along 
£ with large KK masses, the induced energy density is 
exponentially small. In particular, for sufficiently small 
length scales of the internal space this is the case for 
all nonzero KK modes and the main contribution comes 
from the zero mode. For small interbrane distances, to 
the leading order, the induced energy density is given by 
formula l|82|) and the contribution of the induced part 
dominates the single brane parts. When the right brane 
tends to the AdS horizon, 2% — > 00, the induced energy 
density on the left brane vanishes as e~ 2Xl3Zb / z^ 1 ^ 2 for 
the nonzero KK mode and like z b ' Dl ^ 2L ' for the zero 
mode. In the same limit the corresponding energy den- 
sity on the right brane behaves as z ® 2+Dl / 2+1 e -2\,3z b 
for the nonzero KK mode and like z® 2 ~ 2l/ for the zero 
mode. In the limit when the left brane tends to the AdS 
boundary the contribution of a given KK mode into the 
induced energy density vanishes as z^ +2u and as z 2u for 
the left and right branes, respectively. For small values of 
the AdS curvature radius corresponding to strong grav- 
itational fields, for nonzero KK modes under the condi- 
tions \pz a ^> 1 and \(}{zi, — z a ) 1, the contribution 
to the induced energy density is suppressed by the factor 
e -2\ fj (z b ~z a ) _ p Qr zero KK mode, the corresponding 
energy density integrated over the internal space behaves 
as k® 1+1 exp[(Di5j+2v)kD{a — b)] for the brane at y = j 
and is exponentially small. In the model without the in- 
ternal space we have similar behavior with v replaced by 
v\ and for a scalar field with C < Cd+£>i+i the suppres- 
sion is relatively weaker. 

In the model under discussion the hierarchy between 
the fundamental Planck scale and the effective Planck 
scale in the brane universe is generated by the combina- 
tion of redshift and large volume effects. The correspond- 
ing efective Newton's constant on the brane at y = j is 
related to the higher-dimensional fundamental Newton's 
constant by formula l|88|l (see also Eq. I|91[l for the ratio 
of the corresponding Planck masses) and for large in- 
terbrane separations is exponentially small on the brane 
y = b. We show that this mechanism also allows obtain- 
ing a naturally small cosmological constant generated by 
the vacuum quantum fluctuations of a bulk scalar. For 
large interbrane distances the ratio of the induced cos- 
mological constant to the the corresponding Planck scale 
quantity in the brane universe is estimated by formula 
(|93[l and is exponentially small. For the visible brane 
in the higher dimensional generalization of the Randall- 
Sundrum two-brane model, this ratio is given in terms of 
the effective and fundamental Planck masses by Eq. 19611 . 
We have considered two classes of models with the com- 
pactification scale on the visible brane close to the fun- 
damental Planck scale. For the first one the higher di- 
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mensional Planck mass and the AdS inverse radius are 
of the same order and in the second one a separation 
between these scales is assumed. In both cases the cor- 
responding interbrane distances generating the hierarchy 
between the electroweak and Planck scales are smaller 
than those for the model without an internal space and 
the required suppression of the cosmological constant is 
obtained without fine tuning. 

In Section we have considered the total vacuum en- 
ergy in the region between the branes, evaluated as the 
sum of zero-point energies for elementary oscillators. We 
have shown that this energy is equal to the sum of sur- 
face and volume energies. The latter is evaluated as the 
integral of the bulk energy density. Similar to the case 
of the surface energy, for the evaluation of the total vac- 
uum energy we have used the zeta function regularization 
technique. This energy is presented as the sum of the sin- 
gle brane and interference parts. We have extracted the 
pole parts of the single brane vacuum energies and have 
discussed the corresponding renormalization procedure. 
The divergences can be absorbed by adding the countert- 
erms to the brane action which have the structure similar 



to that for the surface energy. The interference part in 
the vacuum energy is given by formula (|127f) and is reg- 
ular for all nonzero interbrane distances. We have inves- 
tigated the asymptotic behavior of this part in various 
limiting cases. For sufficiently small interbrane distances 
the total vacuum energy is dominated by the interference 
part given by formula i|loT|) . In particular, from this for- 
mula it follows that in the case a — 1 any minimum of 
the vacuum energy is local and the fixation of the inter- 
brane distance by the corresponding effective potential 
is stable only locally. Further, we have shown that the 
induced vacuum densities and vacuum effective pressures 
on the branes satisfy the standart energy balance equa- 
tion Q134[l with the inclusion of the surface terms. 
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